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In this paper we propose the recursive stochastic state selection method, an extension 
of the recently developed stochastic state selection method in Monte Carlo calculations 
for quantum spin systems. In this recursive method we use intermediate states to define 
probability functions for stochastic state selections. Then we can diminish variances 
of samplings when we calculate expectation values of the powers of the Hamiltonian. 
In order to show the improvement we perform numerical calculations of the spin-1/2 
anti-ferromagnetic Heisenberg model on the triangular lattice. Examining results on the 
ground state of the 21-site system we confide this method in its effectiveness. We also 
calculate the lowest and the excited energy eigenvalues as well as the static structure 
factor for the 36-site system. The maximum number of basis states kept in a computer 
memory for this system is about 3.6 x 10 7 . Employing a translationally invariant initial 
trial state, we evaluate the lowest energy eigenvalue within 0.5% of the statistical errors. 
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1. Introduction 

Quite recently we have proposed a new method of Monte Carlo calculations, which we call the 
stochastic state selection (SSS) method, ^ to obtain energy eigenvalues of quantum spin systems. 
Under this method we can numerically eliminate a lot of basis states in a mathematically justified 
way. This reduction enables us to calculate expectation values of powers of the Hamiltonian even 
when a limited computer memory is available. In a previous work 2 ) we applied this method to the 
two-dimensional Shastry-Sutherland model and obtained reliable energy eigenvalues for the first 
excited states on a 8 x 8 lattice across the critical region. 

In this paper we propose an extended version of the SSS method, which we will refer to as 
the recursive stochastic state selection (RSSS) method. In the SSS method we have introduced 
random choice matrices, whose elements are stochastic variables to follow probability functions 
named on-off probability functions. Here these probability functions are defined only once at the 
beginning. In contrast to them, the on-off probability functions in the RSSS method are defined in 
sequence for each random choice matrix. By this extension we can expect smaller variances of the 
expectation values of high powers of the Hamiltonian compared to the SSS results. Here we examine 
the effectiveness of the RSSS method in numerical studies of the anti-ferromagnetic Heisenberg 
quantum spin-1/2 system on the triangular lattice. 3 ~ 9 ^ To this typical strongly frustrated system, 
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other methods such as the ordinary Monte Carlo methods or the perturbative calculations are 
hardly applicable. We first show that the method is successful for the model on the lattice of 21 
sites. Then we proceed to the 36-site system to demonstrate that we can estimate the ground state 
energy eigenvalue, a first-excited energy eigenvalue and the static structure factor by this method. 

The plan of this paper is as follows. In the next section we briefly summarize the SSS method 
in order to prepare for the following section and the Appendix. Section 3 is to explain the RSSS 
method, where the recursive procedure of the method is stated. In sections 4 and 5 we apply the 
method to the anti-ferromagnetic Heisenberg spin-1/2 system on the triangular lattice. Section 4 
is devoted to basic examinations. We first examine results on the energy eigenvalue of the ground 
state for a small system of 21 sites. We show that the RSSS results are in good agreement with 
exact values. We also observe much less variances for the expectation values of high powers of 
the Hamiltonian compared to the SSS results. Then we calculate the lowest energy eigenvalue for 
the 36-site system. Here we discuss on a way to find a good approximate state. In addition, we 
state how we estimate the ground state energy eigenvalue from numerically calculated expectation 
values. In section 5 we present further calculations for the 36-site system. Here we employ an 
improved initial trial state which reflects the translational invariance of the Hamiltonian. We 
show that the statistical errors of the ground state energy are much reduced by this improvement. 

0.10 

Our result is E = —20.04 ± , which is about 0.6% higher than the known exact value 

0.06 

E = — 20. 17344. 5 > 7 ' 8 ) We also present results on the energy eigenvalue of an excited state and the 
static structure factor. The final section is devoted to summary and discussions. In the Appendix 
we discuss some properties concerning to the on-off probability functions. 



2. Stochastic State Selection Method 

This section is devoted to a brief summary of the SSS method. ^ Let | ip) be a state of the 
system which is expanded by a basis { | i)} as 

I V>> = £ I i)ci , (l) 
1=1 

Ny being the size of the full vector space. We denote the Hamiltonian of the system by H. In 
order to numerically evaluate the expectation values 

E(L) = (i; | H L | V) (L = l,2,---), (2) 

we define 10 ^ 

E {j]} (L) = (V | HM {T](L)} HM {r)(L - 1)} ■ ■ ■ HM {vW} | V) • (3) 
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The random choice matrices M^( m )j (m = 1, 2, • ■ ■ , L) are the following diagonal matrices 

• \ 

o ,,r' ... o 
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(4) 



with random variables rj^'s. These variables are determined according to the on-off probability 
functions 



max(l,e/|cj|) (q^O) 



e/S 



((H = 0) 



(5) 



where 



P(v,a) = ±6(r)-a)+(l-±) SI,,) 



(6) 



and e and 5(< e) denote given positive constants. Note that each random choice matrix is inde- 
pendently generated using the same on-off probability function Pi(rj) for rf^ m \ We can reduce the 
effective number of the basis by operating M^ m )y to the state ffM^m-ijj - - - HM^i)} I ^) because 
many of generated values of ^"^'s in {r/ m )} = {r/j m ' ) , rj^p , - - - ,T7^} are zero. It is guaranteed 
that the statistical average of E^y(L) gives E(L), 



((E {V} (L))) = E(L) 



(7) 



since the statistical average of -qf^ is 



and {{r}p rip 



■ri\ L) )) = {{rif ) )){{rif ) )) 



POO 

/ »*f m) Pirf m) W m) = 1 , (8) 
{{Vi^)} holds for independently generated {t/ 2 - 1 }, 



3. Recursive Stochastic State Selection Method 

In this section we show a procedure to estimate expectation values following the RSSS method. 
While the on-off probability functions (5) in the SSS method are common to all m for each i, we 
generate the random variables {rj^} in the RSSS method using the preceding intermediate state 
HM {r)(m -i) } ■ ■ ■ HM {r)W} | ip). 



We start from the state | tp^) =| ip), 



N v 



^ (0) ) = Ei^ 0) - 



(9) 



i=i 
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Then we determine the first random choice matrix Mjyi)} generating random variables ryp ac- 
cording to the on-off probability functions 



if >(„>,P („;„?>) . of . J -('.^O , 

oo (c\ ' = 0) 



(10) 



where the function P(rj; a) is given by (6) and e is a positive constant. Calculating the state 

HM {7}W} | V (0) ) with this M {J)W} we obtain | = HM {j){1)} | t/> (0) )/C (1) - Here C^(> 0) is 
given by the equation 

2 



= <^(0) | M {l?(1)} tf 2 M {??(1)} | V {0) ) 



(11) 



which follows from the normalization condition (ip^ \ ip^) = 1. Next we determine the second 
random choice matrix Mr (2)1 in the same manner using 

I V> {1) > = £ I *>Ci 1} instead of | V {0) ), 

namely using the probability functions 



P?\ri) = P(T,;a. 



(2)> 



(2) 
«i = < 



' max(l,e/| Cj {1) |) (cf } / 0) 



(12) 



00 



Then we proceed to calculate ffMujji | C( 2 ) and | ■0( 2 )). Repeating this procedure we 

obtain the m-th normalized state | ip^) from the preceding state | ip( m ~^), 



m) 



= HM {v(m)} \^ m -V)/C< 
C (m) f = (V' (m ~ 1) I M {ri{m)} H 2 M {v ( m)} I ip 1 -" 1 ' 1 ^) . 



(13) 
(14) 



The on-off probability functions to generate M^m+i)} are defined using the coefficients cf 1 ^ in 
I ip (m) ), 



(15) 



as 



,(m+l) 



/ \ _ 7,/ (m+l)\ (m 



+1)= f max(l, e /|cf)|) (cf^O) 



OO 



(cf > = 0) 



(16) 



Thus in the RSSS method we obtain for E(L){ v y defined by (3), 

E(L) {V} = <V (0) I ^C^C^- 1 ) • • • C* 1 ) , 



(17) 



for each set of the random values of {r/} = {{r]^}, ■ ■ ■ , {tj^}}- 

Note that we should be careful here to think of the statistical averages concerning 
because the probability functions P^ (77) 's depend on the precedingly generated sets of random 
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variables {r?W}, {?/ 2 )}, {r^ 1 )}. Yet, we see that {{r]l m) )) = 1 holds even in this case for a 
fixed set of {{^ ( - Tn ~ 1 - ) }, • • • , {v^}}- We therefore obtain the relation (7) if we do recursive samplings 
which ensure 

((((• • • ((((^ (L) ))fixcd{{^-D},...,{r,(i)}} Vj L _1) ))fixed{{^-2)},- ,{„(!)}} • ■■Vk ) ))& X ed{^)} ^f^) = 1 • ( 18 ) 

We expect that variances of ((E^i (L))) in the RSSS method become smaller for large L than 
those in the SSS method because our definition of a^ m+1 ^ in (16) minimizes a quantity S^" 1 ) = 
e 2 ((N^)) + (([g^] )). Here 

g( m ) measures the difference of the truncated state M^ m +i)y | v (m) ) 

from the original one | ip( m ^), 

| x (m) } g (m) = Mwm+1)} | ^M) _ | ^M) ) {x (m) | ^ = X ? (lg) 

and denotes the number of non-zero coefficients in (15). The reason why is minimized 

by a\ m+1 ^ in (16) will be given in the Appendix. 

4. Numerical Examinations 

Here we numerically examine our method. We measure 

^ ^smpl 

«£ W (L))) smpl = E {vh( L ) ' (2°) 

"smpl fc=1 

in order to evaluate E(L), where {rf\k denotes the A;-th set of random values of {{r/ L )}, • • • , {r]^}}. 
In the SSS method, replacing ((E^(L))) by (20) is theoretically allowed because {-q( m ^} is indepen- 
dent of {7/™')} if m / m'. In the RSSS method, on the other hand, things are more complicated 
as was pointed out in the previous section. Since the recursive sampling to ensure (18) is difficult 
to carry out in numerical work because of the bursting sample numbers for large values of L, we 
adopt a sequential generation of only one set of {rf\k = {{r/ L )}, • • • , {r]^}}k for each sampling in 
(20). This replacement is analytically justified in the sense that an iterated integral is replaceable 
by a multiple integral in most cases. We will see our numerical check confirms that we can obtain 
reliable values of E(L) by (20) without recursive samplings. 

After the same manner as the SSS method, we evaluate the error of ((E^(L))) smp i by 



Er(L) 24p^, (21) 

V ^smpl 



where 



p\ v} (L) = « [E {v} (L)] 2 » smpl - [((E {v} (L)» smpl ] 2 



(22) 
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As a concrete example we adopt the anti-ferromagnetic Heisenberg quantum spin system on 
the triangular lattice. The Hamiltonian of the model is 

^ = jE^'<T 3 , (23) 

(m) 

where <Xj is the Pauli matrix on the site i and runs over all nearest neighbor pairs on the 

lattice with the periodic boundary conditions for both directions. 7 ' n ) The coupling J is fixed to 
be 1 throughout this paper. 

First we present in Tables I and II our results for a small lattice with 21 sites. Since we can 
obtain an exact ground state | V'e) by the Lanczos method for this system, we put | i/j^} =| V'e) 
here with = 352,716. Table I shows results on ((E E { v }(L))) smp i up to L = 15 together with 
the exact values Ee(L), where we add a suffix E in order to indicate that an exact eigenstate is 
used as | tp^}. The value of the parameter e is 0.1. In Table II we present ((N® (L))) smp i and 
((N® (£)}}smpi) which denote the numbers of non-zero coefficients before and after operating the 
random choice matrix Mulji to the state HM^l-i)^ ■ ■ ■ HMs^x \ i/ie). For comparison, these 
Tables also contain the results obtained by the SSS method. The value of the parameter e in the 
SSS method is adjusted to be 0.06 so that ((iV^ (L))) smp i nearly equals to that in the RSSS method. 
In fact, Table II shows that ((N^(L))) smpl ~ 0(5 x 10 4 ) for L > 6 in the SSS method, which is 
comparable with ((A^(L))) smpl ~ 0(6 x 10 4 ) for L > 3 in the RSSS method. We see that the RSSS 
results in Table I fairly reproduce the exact values and that, as is expected, errors are less than 
those in the SSS results. Although Er(L) is smaller in the SSS method when L < 3 due to the 
smaller value of e, it becomes smaller in the RSSS method for L > 4. The difference between errors 
in both methods grows as L increases. The RSSS method therefore surely enables us to obtain 
meaningful results for larger values of L compared to the SSS method. 

Now we turn to the larger size system of 36 sites. For our purpose to evaluate the energy 
eigenvalue based on the power method, we first have to find an approximate state for this system. 
Here we use the restructuring techniques 12 ) to form fundamental basis states where three spins 
of each of non-overlapping triangles are diagonalized. Since neither simple way employed in the 
previous work 1 ' 2 ) could not give us a good approximate state, we combine them. Referring to the 
discussions in ref. 13 we start from a trial state and repeat the following procedure until the change 
of the obtained eigenvalue becomes negligibly small. 

( a ) Operate the Hamiltonian to the trial state repeatedly as long as the number of the basis states 
is acceptable to our computer memory resources. Then a truncated vector space is formed. 

( b ) Calculate the lowest energy eigenvalue in the truncated vector space by means of the Lanczos 
method. 

( c ) Drop small coefficients of the obtained eigenstate. Use it as the next trial state. 

Our best approximate state | V'a) obtained from the above procedure is composed of 12,281,253 
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non-zero components. The expectation value of H is (Va I H | Y>a) = — 18.418. 14 ) Using this | tpA) 
we calculate the expectation values. Let Qa{L) = (i/ja \ Q L | V'a), where Q = 51 — H with 
the identity operator /. Since our method is based on the power method, we use Q L instead of 
H L in order to suppress excited states. Results on ((QA{ v y(L))) siav> \ as well as ((-/V^(L))) smp i and 
((iV a A (L))) smpl with e = 0.016 and ra smpi — 10 are presented in Table III up to L — 12. Here 
the maximum value of L is determined under a condition that the square root of the variance, 



p^y(L), is less than the expectation value. The maximum program size is about one Gbytes 
and the CPU time is about one hour for one sampling by a Pentium 4 machine. The results in 
Table III show that the RSSS method is effective, not only for | Ve) of the 21-site system but also 
for | ^a) of the 36-site one, to suppress a rapid increase of the variance for large L. For example, 
the variance pf^ity m the SSS method is 9.3 x 10 13 , which is about 600 times as large as the 
RSSS result P{ v y{5) = 1-6 x 10 11 . We also see in Table III that the number of non-zero coefficients 
drastically reduces by the random choice matrix, ((N^(L))) smp \ <C ((N£-(L))) smp \ for each L, and it 
increases very slowly as L grows. In Figure 1 we plot the ratio ((QA{r/}(L))}smp\/((QA{ri}(L — l))) S mpi 
(((£?A{»7}(0)))smpi = 1), which should give the exact value 5 — E in the large L limit. We observe 
that the ratio increases for the data up to L = 10. This increase suggests that the power method 
works. For L = 11 and 12 the errors are too large to see this tendency. We therefore exclude 
((QAMCU^smpl and ((<9 A {r,}( 12 )))smpl in the following analysis. 

In order to estimate the energy eigenvalue from these expectation values we use the same fitting 
form as was introduced in the previous work, 1 ) 



W»a| Wa> = Q L Uo + * ' = F(L, Q, q , gi, a) , (24) 
V L + a + 1/ 

go, Qi, ol being free parameters which are to be determined together with Q by the fit. Changing 
these four parameters under two constraints go+9i/( a +l) = (^A I V'a) = 1 and Q{go+gi/ (a+2)} = 
(VaJQIVa.) = —18.418, we look for the minimum of the difference D, 

n-V'L ((QA{r,}(L))) smp i ] 2 

tA F(L,Q,g , qi ,a)\ ' [2b) 
with L max = 10. We accept the values of the parameters if D is less than the sum of the relative 



errors 



E 

L=2 



Er(L) 



1 2 



((QA{ V }(L))) S m P \_ 



2.2 x 10~ 3 



1.37 

The result of this fit is E& t = 5 — Qm = —20.50 ± . The function F(L, Q& t , </ofit 5 Qiftti a &t) 

0.12 

with fitted parameters to minimize D is plotted with a solid line in Fig. 1. 
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5. Further Calculations 

In this section we add several results for the 36-site system. In the previous section we obtained, 
starting from an approximate state | V'a}) the ground state energy of the 36-site system which is in 
agreement with the known exact value. Statistical errors, however, are not satisfactorily small. One 
way to decrease the errors is to increase the number of samples, but it is quite time-consuming. 
In this section we employ another way, which is to improve the trial state by requesting the 
translational invariance. For this purpose we introduce a wave vector k. Roughly every N s basis 
states are linearly combined to form one new basis state for a definite value of k, where N s denotes 
the number of sites of the system. Requesting k = we construct an improved trial state | i^t), 
whose number of new basis states amounts 13,911,394 and (t/>t I H | t/>t) = —19.710. Here we 
abandon the restructuring techniques in order to avoid too much complexity. Results up to L = 6 
from 140 samples with e = 0.02 are shown in Table IV and Fig. 1. Using the same fitting form (24) 

we obtain = —20.04 ± , whose statistical errors are less than 0.5%. The number of non- 

0.06 

zero coefficients before (after) operating the 6-th random choice matrix is ((-/V^(6))) smp i ~ 3.0 x 10 7 
(«iV a T (6)» smpl ~ 7.9xl0 5 ). 

The lowest energy eigenvalues for excited states with fixed values of S z , the z-component 
of the total spin, are also calculable in similar manners. Our results for the S z = 1 homoge- 
neous (k = 0) state are presented in Table V. Here e = 0.02 and n smp i = 88 for L < 5, while 
e = 0.03 and n smp i = 180 for L > 6. We start from a trial state | 4>t{S z = 1)), for which 
(4>t(S z = 1) | H | 4>t{S z = 1)) = —18.813. Up to L = 5 we, as usual, calculate inner products 
of (iPt(S z = 1) | and QM^ L )y ■ ■ QM^i)^ \ ipT(S z = 1)). For L > 6, on the other hand, we 
measure inner products between independently generated ('4>t(S z = 1) | Mjy^^Q • • • M^^^^^Q 
and QM^ v (5)y ■ ■ ■ QM^i)^ \ iPt(S z = 1)) in order to save memory resources. This measurement is 
allowed because (( (V | M {r > W} Q ■ ■ ■ M {r] , (L _ 5)} Q ■ QM {t;(5)} • ■•Q%i )} | V) )) = (( <V> I Q L I 4>) )) 

1*1 / °- 27 

holds. 10 - 1 Based on the assumption (24), we obtain E^ t = —19.21 ± from data in Table V. 

0.10 

Finally, in order to demonstrate that we can extract physical properties of the ground state 
with this method, we report results on the static structure factor F(k) defined by 

F(k) <V> E | S_ k • S k | ^ E ) , S k ee — L= ]T aje-^ , (26) 



2v^ 

where j runs over all site of the lattice and rj denotes the position vector for site j. In numerical 
study we define 

ftfr | Q m s_ k • s k Cr I til (27) 
(V't I Q 2m I V-t) 
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which gives F(k) in the large m limit, 

lim F m (k) = F(k) . 

m^oo 

Then we evaluate F m (k) from measurements of 

(( (V't | QM {v(2m -2) } ■ ■ ■ QM {v(m)} S_ k • S k QM {v(m - 1)} ■ ■ ■ QM {r)W} I Vt) })sm P i 
(( (V't | QM {r p m - 2) y ■ ■ QM {v(m)} QM {v(m -i) } ■ ■ ■ QM {vW} \ ip T ) )) smp i 

Our results for e = 0.03, n smp i = 8 and m = 6 are shown in Fig. 2. Here we do not need many 
samples because we allow maximally one per cent of the statistical errors for F m (k). They are 
in good agreement with the data from ref. 7 except for those with k/ko = 1. When k = ko we 
can obtain a better value from the squared sublattice magnetization Ai 2 , using the relation 16 ) 
2A r s F(ko) = 9.M 2 . The reason for this better value is that in measurement of M 2 one does not 
suffer from heavy cancellations. We find this value is much closer to the datum given in ref. 7. 

6. Summary and Discussions 

In this paper we developed the recursive stochastic state selection (RSSS) method, which is 
derived from the stochastic state selection (SSS) method. 1 ) As the word "recursive" indicates, 
on-off probability functions to generate the random choice matrix are recurrently determined so 
that they reflect the newest intermediate state at each step of the RSSS procedure. The merit of 
this modification is that variances for expectation values with high powers of the Hamiltonian grow 
slowly compared with those in the SSS method. 

As a result of less variances we can now evaluate the ground state energy for the triangular 
lattice Heisenberg spin- 1/2 system by means of the stochastic selections. On a relatively small 
21-site lattice we examined the RSSS method using the exact eigenstate. Results there show that 
the variance p^(10) was reduced to ~ 1/2000 of the one in the SSS method. On a larger 36-site 
lattice we calculated the expectation values of (5/ — H) L from two approximate states | V'a) an d 
| ^t), where the latter realizes the wave vector k = 0. We again emphasize that the numerical 
study of such high powers would not become possible without introducing the RSSS method. From 
obtained data with | V"r) we estimated the energy of the ground state E based on the extrapolation 
assumption. The result indicates —20.10 < E < —19.94, which is slightly higher than the known 
exact value —20.17344. On the 36-site triangular lattice we also calculated expectation values on the 
lowest energy eigenvalue with the S z = 1 state. Further, we calculated the static structure factor, 
which is a typical measurement to obtain physical properties of the ground state. Our results fully 
ensure that not only energy eigenvalues but also the static structure factor is calculable in our 
method. 

Let us finally comment on whether it is possible to study larger systems, the 48-site system 
for instance, by this method. In a simple estimation by a guess that the cost is proportional to 



(28) 
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the system size, we would need 2 48 /2 36 ~ 4 x 10 3 times of computer resources compared to the 
the CPU time and the memory we used in our present study for the 36-site system. This factor is 
too large for us to easily overcome. We, however, show in this study that a better trial state with 
a fixed k can much reduce the cost of calculations. We hope, therefore, we can make a numerical 
study of larger triangular spin systems in the near future. 

Appendix: On-Off Probability Functions 

In this Appendix we consider a random choice matrix Mr^i = diag.{r/i, 772 , • • • ,tjn v }, where 
each random variable 77, is generated by an on-off probability function defined with a parameter £j 



Our purpose is to discuss what values of {£} = {£1, £2, • • • , 6v v } are favorable in order to decrease 
the variances of the expectation values. We will see that the definition of the on-off probability 
functions in the RSSS method is "best" in the sense shown below. 

First, let Ni denote the number of non-zero r]i, which is therefore 1 (0) when r/j / (r/j = 0). 
Following the discussion in ref. 1, we obtain for the statistical average of N = ^ Ni in the above 





(A-l) 





(A-2) 



because 




(A-3) 



P Ni (Ni) = T S (Ni ~ 1) + (1 - hs(Ni) 




(A-4) 



Next, multiplying M^y to a normalized state vector | \I>), 




(A-5) 



we obtain 





(A-6) 




(A-7) 
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where the normalization condition | VP) = ^ c 2 = 1 is used. Let us then introduce a state vector 
| x)g which represents the difference between | vp) and | VP), 

| X )9 = M {v} | *)- | *) = I i)<H (Vi - 1) , (A-8) 

% 

with a normalized | x). It should be emphasized here that (( (<£ | x)9 )) = holds for any state 
| = X) I because ({r]i)) = 1. From (A-8) we see 

(( (* | [M W ] 2 | *) )) = (( (vi/ | vi/) + 2g( X | *) + g 2 ( X \ x ) )) = 1 + « 5 2 )). (A-9) 

Using 

/>oo 

((^ 2 » = / ri 1 i P i { m )d m = (, i (A-10) 

JO 

we obtain from (A-7) and (A-9), 

((g 2 )) = «E c * 2 W -!))) = E c ? «» - 1) = E c ? (6 - 1) • ( A - n ) 

i i i 

Is it possible to find any {£} which makes both of {(N)) and ((g 2 )} small? From (A-2) and (A-11) 
we learn that we cannot lessen ((g 2 )) without increasing ((N)). Let us consider, then, to minimize 
a quantity 

AT V , 2 . 

e 2 ((N)) + ((g 2 )) ^j-^te-Dj^^,..,^), (A-12) 

with a positive constant e in ranges £j > 1 for all i. Since dS/d^i = cf — e 2 /£ 2 , we find 

• dS/d^i = at £ = e/|cj| (> 1) if < |q| < e. 

• dS/d^i > for £j > 1 when |cj| > e. In this case £j = 1 gives the minimum of e 2 /£j + c 2 (£j — 1). 

• dS/d^i < if Ci = 0. In this case e 2 /& + c 2 (& - 1) = e 2 /& -► when & -► oo. 
Therefore the minimum value of S 1 £2, ■ ■ ■ 5 6v v ) is realized by 

^ = I max(l,e/|ci|) (q / 0) g 
[00 (q = 0) 

To summarize this Appendix, we find the "best" choice of {£} in (A-l), by which the quantity 
e 2 ((Af)) + ((g 2 )} in (A-12) is minimized. 
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15) For example, using (19) we see that 

QM { ^ )} \i>) = Q\4>) + Q\x)g, 

Q M { V '(i)} I V>> = Q\i>) + Q\x'W , 

and then 

« (i,\M W(1)} Q-QM {vll)} \i>} )) 

= (( <v I Q 2 ! V) )) + « g'ix \Q 2 \^))) + « <^ ! Q 2 1 x)g )) + « g'ix I Q 2 1 x)g )) 
= «<^|Q 2 j^))). 



16) This is an identity for systems with the total spin S = 0. 
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RSSS 




SSS 




L 


Ee{L) 






(( £: E{r ? }(^)))smpl 




1 


-0.117809 


xlO 2 


- (0.117787 ± 0.000072) 


xlO 2 


-(0.117796 ± 0.000050) 


xlO 2 


2 


0.138790 


xlO 3 


(0.13885 ± 0.00015) 


xlO 3 


(0.138760 ± 0.000098) 


xlO 3 


3 


- 0.163507 


xlO 4 


-(0.16362± 0.00027) 


xlO 4 


-(0.16343 ± 0.00020) 


xlO 4 


4 


0.192626 


xlO 5 


(0.19286 ± 0.00044) 


xlO 5 


(0.19244 ± 0.00067) 


xlO 5 


5 


-0.226931 


xlO 6 


-(0.22710 ± 0.00070) 


xlO 6 


-(0.2270± 0.0030) 


xlO 6 


6 


0.267346 


xlO 7 


(0.2671 ± 0.0011) 


xlO 7 


(0.269± 0.015) 


xlO 7 


7 


- 0.314958 


xlO 8 


-(0.3152 ± 0.0018) 


xlO 8 


-(0.286± 0.049) 


xlO 8 


8 


0.371049 


xlO 9 


(0.3708 ± 0.0030) 


xlO 9 


(0.24± 0.17) 


xlO 9 


9 


- 0.437129 


xlO 10 


-(0.4357 ± 0.0050) 


xlO 10 


- (0.39± 0.17) 


xlO 10 


10 


0.514978 


xlO 11 


(0.5170 ± 0.0085) 


xlO 11 


(0.30± 0.41) 


xlO 11 


11 


- 0.606691 


xlO 12 


- (0.607 ± 0.014) 


xlO 12 


— 




12 


0.714737 


xlO 13 


(0.709 ± 0.025) 


x 10 13 






13 


- 0.842025 


xlO 14 


-(0.850 ± 0.042) 


xlO 14 






14 


0.991982 


xlO 15 


(1.005 ± 0.072) 


xlO 15 






15 


- 0.116865 


xlO 17 


-(0.122 ± 0.013) 


xlO 17 







Table I. Results on ((E E { v y(L))) Blnp \ (L = 1, 2, • • • , 15) obtained for the 21-site triangular lattice by the RSSS method 
with e = 0.1 and the SSS method with e = 0.06. The number of samples is 10 4 . Exact values of Ee(L) are also 
presented for comparison. 
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RSSS 


sss 


L 


((nP(l))) , 


((N E (L))) i 

\\ iv a ^-^//smpl 


((N*(L))) smpl «iV a E (L)» smpl 


1 


352716 D 


2659 3 


352716.0 4432.2 


2 


74254.9 


2038 3 


108023.1 2603.0 


3 


60343 5 


1923.6 


70882.0 2197.7 


4 


57531.4 


1920.4 


60150.5 2048.7 


5 


57365 5 


1946.4 


56182.7 1984.9 


6 


57819 9 


1975.7 


54473.2 1955.8 


7 


58330 5 


2000 4 


53687.7 1941.7 


8 


58775 9 


2020 2 


53304.5 1934.3 


9 


59672 5 


2033 5 


53102.0 1930.7 


10 


60088 4 


2043.6 


53003.4 1929.1 


11 


60185 8 


2050.4 




12 


60143.4 


2053 9 




13 


60179.8 


2057.8 




14 


60624.2 


2059.9 




15 


60691.2 


2061.4 





Table II. Numbers of non-zero coefficients before (((A^(L))) smp i) and after (((N® (L)}) smp i) operating the random 
choice matrix M^Dy to the state HM^l-i) y ■ ■ ■ HM^ijy \ ips) for the 21-site triangular lattice obtained from 
10 4 samples with e = 0.1 (RSSS) or e = 0.06 (SSS). 
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L 


((QA{ v }( L )))smpl 




((K(L))) smpl 


«^)»smpl 


1 


(0.234137 ± 0.000077) 


xlO 2 


12281253.0 


91378.8 


2 


(0.55286 ± 0.00045) 


xlO 3 


21738077.0 


69452.6 


3 


(0.13114 ± 0.00019) 


xlO 5 


19123511.3 


69429.2 


4 


(0.31260 ± 0.00071) 


xlO 6 


19821436.4 


74755.3 


5 


(0.7469 ± 0.0025) 


xlO 7 


21556264.4 


81933.1 


6 


(0.17941 ± 0.00093) 


xlO 9 


23674697.9 


89782.3 


7 


(0.4296 ± 0.0035) 


xlO 10 


25924492.7 


97767.5 


8 


(0.1028 ± 0.0014) 


xlO 12 


28180669.0 


105616.8 


9 


(0.2507 ± 0.0055) 


xlO 13 


30379138.4 


113156.1 


10 


(0.629 ± 0.024) 


xlO 14 


32478107.9 


120316.8 


11 


(0.153 ± 0.010) 


xlO 16 


34462103.9 


127006.7 


12 


(0.355 ± 0.042) 


xlO 17 


36309064.0 


133239.7 



Table III. Results on ((Q A {„}(£)»smpi, ((A^(L))> smp i and «iV a A (I0)) smp i (L = 1, 2, • • • , 12), obtained for the 36-sitc 
restructured triangular lattice by the RSSS method with e = 0.016, where Qa(L) = (i/>a | {51 — H) L \ ipA) with 
the identity operator I. The number of samples is 10 3 . 
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L 


((<2T{r,}(£)))smpl 




«iV b T (L)» smpl 


«iV a T (L))) smpl 


1 


(0.247073 ± 0.000045) 


xlO 2 


13911394.0 


398861.8 


2 


(0.612899 ± 0.000193) 


xlO 3 


13513630.8 


437278.8 


3 


(0.152364 ± 0.000063) 


xlO 5 


16928438.7 


518694.1 


4 


(0.37929 ± 0.00018) 


xlO 6 


21395583.2 


610462.3 


5 


(0.94546 ± 0.00057) 


xlO 7 


25728449.4 


701952.5 


6 


(0.23594 ± 0.00017) 


xlO 9 


30450059.3 


787238.1 



Table IV. Results on ((Q T{v} {L))) smph ((Ng (L)» S m P i and «iV a T (L))) smpl (L = 1, 2, • • • , 12) obtained for the 36-site 
triangular lattice by the RSSS method with e = 0.02, where a trial state | Vt) is used instead of | iPa)- The number 
of samples is 140. 



28 



26 



24 

CC 

22 
20 



18 

5 10 15 

L 

Fig. 1. Ratios R A = {{Q A { v }(L)}} S m P i/{{QA{ v }{L- l)» S m P i (crosses) and R T = ((QT{ v }(L)}} smp i/ {{Qt{ v } (L - l))) smp i 
(pluses) for the Heisenberg spin model on the 36-site triangular lattice. Values of Ra are obtained from the results in 
Table III and {{Qam (0))) smp i = 1- Values of _Rt are obtained from the results in Table IV and {{Qtm (0))) S m P i = 1. 
The dashed line indicates the exact eigenvalue of Q = 5 — _E. 5 ' 7 ' 8 ' The solid lines present fitted values of (24). 
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L 


((QT{ V }(L))}smpl 




«^b T W»smpl 


«iV a T (L)» smpl 


1 


(0.238085 ± 0.000067) 


xlO 2 


13957843.0 


557850.5 


2 


(0.569663 ± 0.000272) 


xlO 3 


17375161.1 


553892.4 


3 


(0.136574 ± 0.000096) 


xlO 5 


20943597.7 


618879.1 


4 


(0.32779 ± 0.00026) 


xlO 6 


25154518.3 


699008.8 


5 


(0.78765 ± 0.00079) 


xlO 7 


29199658.1 


779763.9 


6 


(0.18980 ± 0.00025) 


xlO 9 






7 


(0.45699 ± 0.00064) 


xlO 10 






8 


(0.11010 ± 0.00017) 


xlO 12 






9 


(0.26527 ± 0.00049) 


xlO 13 






10 


(0.6401 ± 0.0015) 


xlO 14 







Table V. Results on ((Q T{l)} (L))) smp i for L = 1, 2, • • • , 10 and ((Ng (L)» Bmp i and ((7Vj(L))) Bmp i for (L = 1, 2, • • • , 5) 
in the S z = 1 case, obtained for the 36-site triangular lattice from a translational invariant trial state | iPt{S z = 1)}. 
Parameters are e = 0.02 and n smp i = 88 up to L = 5, and e = 0.03 and n amp i = 180 for L > 6. 
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N 




0.2 0.4 0.6 0.8 1 

k/kO 



Fig. 2. The static structure factor Z = 8F(k)/(A r s + 6) for the 36-site system versus k/ko, where ±ko are the two 
wave vectors of the corners of the crystal Brillouin zone, k = |k| and N s (= 36) is the number of sites. Asterisks 
combined by lines to guide eyes are data from ref. 7. Our results obtained by 8 samples of F6(k) with e = 0.03 are 
plotted by open diamonds. A filled square is the value calculated by F(ko) = 4.5M 2 /N B , where M 2 denotes the 
squared sublattice magnetization. 



